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Abstract. Comprehensive VLBI and multi-waveband monitoring indicate that a single superluminal knot can cause 
a number of 7-ray flares at different locations. However, the often very rapid variability timescale is a challenge to 
theoretical models when a given flare (perhaps the majority of those observed) is inferred from observations to lie 
near the 43 GHz core, parsecs from the central engine. We present some relevant observational results, using the 
BL Lac object AO 0235+164 as an example. We propose a turbulent cell model leading to a frequency-dependent 
flUing factor of the emission region. This feature of the model can provide a solution to the timescale dilemma and 
other characteristics of blazar emission. 



1. Introduction 

There exists a crisis in the interpretation the results of 
multi-waveband monitoring of blazars. Comparisons of 
dates of peak 7-ray flux measured by EGRET with epochs 
of ejections of superluminal radio knots (Jorstad et al. 
I2001ap . onsets of millimeter- wave flares (Lahteenmaki & 
Valtaoia l2003p . and changes in centimeter- wave polariza- 
tion (Jorstad et al. I2001b|) lead to the conclusion that 
most 7-rays outbursts are coincident with radio events. 
Yet we know that the jets of blazars are opaque to radio 
emission within ^ 1 pc or more from the central engine. 
This is difficult to reconcile with the shortest time scales of 
GeV variability observed — a few hours in PKS 1622—297 
(Mattox et aL [TMf|) and 3C 454.3 (Foschini et aL l^OTU)) . 
However time scales of variability reflect size scales, 
not distance from the central engine. Furthermore, jets are 
very narrow, with opening half-angle ~ 10°/r, where F is 
the bulk Lorentz factor of the flow in the jet (Jorstad et 
al. I2005p . Because of this, a jet can have a cross-sectional 



radius much less than its distance from the central engine. 
Nevertheless, this does not solve the entire problem. The 
core on 43 GHz Very Long Baseline Array (VLBA) images 
of PKS 1510—089 has been estimated to lie ~ 20 pc from 
the central engine (Marscher et al. I2010a| . Offsets of par- 
secs are implied for other luminous blazars as well (e.g., 
Marscher et al. [MI51 Chatterjee et al. [M?5)) . The cross- 
sectional radius of a luminous blazar is ~ 0.1 pc near the 
43 GHz core. With a Doppler factor (5 ~ 20 and redshift 
~ 0.5, the shortest time scale of variability should be ^ 1 
week, not a fraction of a day. 

Because of this discrepancy, a number of authors have 
insisted that the 7-ray emitting region lies within ~ 10^^ 
cm of the central engine (e.g., Tavecchio et al. l2010p . This 
has the advantage for high-energy emission models that 
optical-uv photons from the broad emission-line region are 
available for scattering to 7-ray energies by highly rela- 
tivistic electrons in the jet. Poutanen & Stern (120101) find 
that pair production off the uv photon field can then cause 



the observed sharp break in the spectrum of some blazars 
at GeV energies. However, such an interpretation requires 
that the timing of the radio and 7-ray events is a chance 
coincidence in every case. This seems unlikely, especially 
with singular events such as the ultra-high-amplitude op- 
tical flare, very sharp 7-ray flare, and passage of a superlu- 
minal knot through the 43 GHz core in PKS 1510-089 at 
essentially the same time (Marscher et al. I2010all2010b|) . 
Furthermore, the number of 7-ray outbursts observed with 
Fermi that are either coincident with or follow the pas- 
sage of a new superluminal knot through the centroid of 
the core is becoming high enough to conclude that a large 
fraction of the 7-ray emission originates parsecs away from 
the central engine (Jorstad et al. 120101 and these proceed- 
ings). 

We propose to solve this dilemma by developing a 
model in which much of the optical and high-energy radi- 
ation in a blazar is emitted near the 43 GHz core in VLBA 
images, parsecs from the central engine. The model allows 
for short time scales of optical and 7-ray variability by 
restricting the highest-energy electrons radiating at these 
frequencies to small sub-regions of the jet. That is, in our 
model the filling factor at high frequencies is relatively 
low, while that of the electrons radiating at ~ 10^°"^'^ Hz 
is near unity. Such a model is consistent with other promi- 
nent features of optical vs. lower frequency emission. 

2. Properties of blazar emission to incorporate in 
the development of a model 

We now have at our avail light curves at many frequencies 
across the electromagnetic spectrum, polarization vs. time 
data at radio, mm-wave, and optical/near-IR bands, and 
well-sampled sequences of VLBI images — in both total 
and polarized intensity — at high radio frequencies for 
many blazars. This is providing a wealth of information 
that must be considered when developing models for the 
nonthermal emission from blazar jets. This information 
includes the following: 
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Fig. 1. Multi- waveband light curves of the BL Lac object 
AO 0235+164 in late 2008. The time when the superlu- 
minal knot (see Fig. 2) passed through the core is marked 
by a vertical line, with the horizontal line in panel 3 indi- 
cating the uncertainty in the date of this event. 



(1) In most blazars, the detailed properties of the 
synchrotron emission change gradually as the wavelength 
of observation decreases from the millimeter to optical 
range, (a) The spectral index (a = —d logF^/d logv) 
steepens from a value usually in the range of 0.5-1.0 
to > 1; (b) the linear polarization increases on average 
and becomes more highly variable (D'Arcangelo et 
al. I2009P : and (c) the time scale of flux variability, 
[ivar = (^2 — ii) I \ln{F2 / Fi)\, where ti and t2 are times 
when flux densities Fi and F2 are measured, respectively] 
decreases. An example of the latter is the BL Lac object 
AO 0235-1-164 [z = 0.94) during a multi-waveband 
outburst in 2008 (Fig. [1]). As seen in Figure [21 a bright, 
extremely superluminal (^ 70c!) knot passed through the 
core at essentially the same time as the peak of a sharp 
optical flare during a ^ 80-day 7-ray outburst. During 
the declining phase of another major flare observed at 
two optical and two near-IR bands, the shortest measured 
time scale of flux decrease depended on frequency u as 



U 



ex V 



-0.16 



^see Fig. [3]). 



(2) The period of elevated optical flux in 0235-1-164 
(see Fig. [1|) matches that of the high 7-ray state, but 
there were no high-amplitude 7-ray flares during the two 
main optical events. Yet electrons with sufficient energies 
to emit optical synchrotron radiation for typical magnetic 
field strengths of - 0.5 G (Hagen-Thorn et al. [2008)) . 
E ~ W^mc^, should also be able to scatter those photons 
(or other optical photons from the surroundings) to 7-ray 
energies. 

(3) The power spectral density (PSD) of flux variations 
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Fig. 2. VLBA images at 43 GHz of 0235+164 in late 2008. 
Contours (factors of 2 starting at 0.023 Jy/beam) corre- 
spond to total intensity and gray scale to polarized inten- 
sity. Line segments indicate direction of polarization elec- 
tric vectors. A circular Gaussian restoring beam is used 
with FWHM similar to the resolution of the longest VLBA 
baselines. Motion of the superluminal knot is indicated by 
the diagonal line. 

of the 7-ray bright blazar 3C 279 follows a power law 
(Chatterjee et al. 12008 ) without any indication of a break 
on time scales from ~ 1 day to 12 years. 

(4) The 43 GHz core on VLBI images lies at a sig- 
nificant distance from the central engine. This must 
be the case theoretically based on the self- absorption 
turnover in the synchrotron Fi, vs. u spectrum at mil- 
limeter wavelengths (Impey & Neugebauer I1988p . If 
emission in the jet continued in a self-similar way down 
to a size less than about 10^^ cm, the turnover would 
occur at much shorter wavelengths. In the radio galaxy 
3G 120, time delays between a drop in X-ray flux from 
the central engine and the passage of a new superluminal 
knot through the core at 43 GHz led Marscher et al. 
([2002)) and Chatterjee et al. ([2009)) to conclude that the 
core lies 0.5 pc from the black hole. In a more powerful 
source, such as the quasars 3C 279 and PKS 1510—089, 
timing studies indicate that the core lies tens of parsecs 
(after de-projection) from the central engine (Marscher 
et al. I2010al [20T0bl Chatterjee et al [2008)) . 

3. A turbulent cell model for variability in blazars 

The red-noise behavior of the PSD of 3C 279 (Chatterjee 
et al. I2008P implies that the variations are stochastic in 
nature. If main events in the light curves associated with 
bright superluminal knots are caused by shocks, these 
shocks must be produced randomly with a range of com- 
pression ratios. The short-term fluctuations can be under- 
stood as the consequence of turbulent ambient jet plasma 
that passes through the shocks (Marscher, Gear, & Travis 
I1992p . If this is the case, then the multi- wavelength be- 
havior of the emission depends on the range of physical 
parameters of the turbulence. 

We propose a model in which the emission arises from 
a collection of turbulent cells of similar size, each with 
strength of the magnetic fleld B (with random orienta- 
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Fig. 3. Top: Multi-band optical/near-IR light curves 
of 0235-fl64 during the interval displayed in Fig. [T] 
Magnitudes here are not corrected for extinction. Bottom: 
Time scale of flux variability (as defined in the text) as 
a function of optical/near-IR frequency of 0235-1-164 be- 
tween Juhan dates 2454756.7 and 2454758.7. Diagonal 
line represents a power law with slope —0.16. Data 
are from the Yale Fermi/SMARTS project, website 
|http://www. astro. yale.edu/smarts/glast/. 



tion) and density of electrons ric that follow a probabil- 
ity distribution, with higher values less likely. We imagine 
that the turbulent cells pass into and out of a disk-shaped 
primary emission region representing a shock, bounded by 
the shock front and, probably, a rarefaction. Electrons are 
injected with a power-law energy distribution at the shock 
front, after which the cell advects away from the shock. 
The cooling time of the electrons depends on B. If inverse 
Compton losses are important, then it will also depend on 
the energy density of photons n-y] we defer consideration 
of this case to a subsequent study. 

The key addition to this model over that proposed in 
Marscher et al. (|1992p is that the maximum electron en- 
ergy i?max varies from cell to cell. According to models of 
energization of particles at relativistic shock fronts (e.g.. 
Baring I 2010p . the efficiency of particle acceleration de- 
pends inversely on the ratio of the gyro-radius to mean 
free path for pitch-angle scattering. This ratio should be 
larger when B is stronger, so that the amplitude of the 
magnetic turbulence rises. If the value of i? in a cell follows 




Fig. 4. Sketch of our proposed model consisting of a num- 
ber of emission cells (circles) inside a disk (represent- 
ing, for example, a shock wave), each with a different 
maximum energy of relativistic electrons. All numbered 
cells have electrons with sufficiently high energies to emit 
synchrotron radiation at a relatively low frequency, e.g., 
< 10^^ Hz. Only the cells marked 2 and higher have elec- 
tron energies high enough to emit at > 10^^'^ Hz. Only 
those marked 3, 4, and 5 can emit at > lO""^^ Hz, only those 
numbered 4 and 5 can emit at > 10^"''^ Hz, and only those 
numbered 5 can emit at > 10^^ Hz. 



a power-law probability distribution above some minimum 
value, then we might expect i?max to do the same. (Or the 
slope of the resultant electron energy distribution might 
be flatter for stronger fields, thereby increasing the relative 
number of electrons with the highest energies. This would 
produce a similar effect on the synchrotron spectrum as 
the case considered here.) 

The different values of B, -Emax, and tic in differ- 
ent cells, with higher values occurring more rarely, leads 
to a frequency dependence of the synchrotron emission. 
Recall that, in the scenario envisioned by Marscher & Gear 
(|1985p . the synchrotron spectrum steepens by 1/2 because 
the dimension transverse to the shock front, and therefore 
the volume V, of the emitting region follows V{iy) ex v~^/'^. 
This will still apply at sufficiently high frequencies if elec- 
trons are injected only as a cell passes across the shock 
front. However, we introduce an additional frequency de- 
pendence on the volume because only some fraction of the 
cells contain electrons that can radiate at frequency v af- 
ter they cross the shock front. We illustrate the geometry 
of our model in Figure SI 

We have not yet developed a definite relationship be- 
tween the electron energy distribution and the power spec- 
trum of the turbulence. Our model therefore does not 
at this stage predict the slope of the spectrum at fre- 
quencies where V^v) < 1. We can, however, turn to ob- 
servations to infer the required functional form of V{h'). 
The relationship tvar oc i/~'^-^^ found above for 0235-1-164 
could correspond to stochastic variations in N cells if 
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Fig. 5. Optical/near-IR colors vs. B magnitude of 
0235+164 during the interval covered in Fig. [T] All mag- 
nitudes and colors have been corrected for foreground red- 
dening in both our Galaxy and an intervening galaxy at 
redshift 0.524 according to Table 5 of 2005, Comparison 
with the horizontal dotted lines shows the trend of flat- 
ter spectra as the flux increases. For example, the B — J 
color difference between 2.4 dX B — 14.17 and 2.85 at 
B — 2.85 corresponds to a change in the spectral index 
from a = 1.01 to 1.34. Data source is the same as for Fig. 

m 



N{v) oc V{v) oc zy^"'^^, since we expect that the level 
of fluctuations in flux (AF/F) oc iV~^/^. This is in agree- 
ment with the steepening of the continuum spectrum from 
far-IR {a w 1.0) to optical {{a « 1.3) at times when there 
is no flare. Since 0.32 is less than 0.5, we infer that cells 
can cross the entire emission region before electrons that 
radiate at B band suffer significant synchrotron losses in 
0235-M64. 

Since a rapid flare corresponds to one or more new cells 
with higher than average B or rie, the spectrum should 
flatten near the peak of the flare, since the volume of the 
enhanced emission is roughly the same at all wavelengths 
within a relatively narrow range (e.g., optical/near- IR) . 
As Figure \5\ demonstrates, this indeed occurred during 
the 80-day 7-ray outburst of 0235-1-164. The change in 
spectral index closely matches that required by the model 
(Aa = 0.32). 

The higher mean level and amplitude of fluctuations 
of linear polarization at higher frequencies is a natural 
consequence of the turbulent cell model. In the extreme 
case of either very little compression at the electron injec- 
tion front (weak shock) or a shock viewed almost exactly 
face-on (so that the component of B in the plane of the 
sky remains chaotic), the mean degree of polarization at 
frequency v would be ~ 0.75N{i>)~^/^ (Burn I966j) and 
the standard deviation - 0.75[2Niu)]-^/'^ (Jones [HHSD- 
This becomes more complicated when the compression is 
substantial and the front lies at an angle to the sky plane. 
There will still be frequency- dependent polarization, but 
simulations are needed for a more quantitative analysis. 

The time scale of variability ivar c>c v~'-''^^ in 0235-1-164 
at optical/near-IR frequencies, probably extending down 
to ~ lO^'^ Hz. If we assume that the volume of emission 
at this frequency includes the entire superluminal knot, 



whose angular diameter was 0.05 mas on JD 2454762 as 
measured by modeling the feature as a face-on disk and 
fitting the VLBA data. This corresponds to a radius of 0.2 
pc. If we adopt 6 « 70, the typical time scale of variability 
at 1/ ^ 10^'^ Hz should be ~ 7 days. This becomes ~ 3 days 
at B band, close to the observed value of 4 days (see Fig. 

El). 

Inverse Compton emission depends on ric and n-y, while 
the synchrotron flux depends on the former and on the 
strength and direction of the local magnetic field. We 
therefore expect correlations between the optical and 7- 
ray emission, but not exact correspondence of light curves. 
The geometrical break in the synchrotron spectrum should 
appear also in the 7-ray spectrum from inverse Compton 
scattering. We plan to investigate this possibility alongside 
the detailed time behavior of the multi- waveband emission 
through numerical simulations of the model in the near fu- 
ture. 
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